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ABSTRACT 

By combining data from different published 3-D simulations of Keplerian 
shearing boxes unstable to the magnetorotational instability (MRI), we highlight 
tight anti-correlations between the total effective inferred angular momentum 
transport parameter, a to t, its separate Maxwell and Reynolds contributions a mag 
and akin, and the kinetic to magnetic pressure ratio /3, defined with the initial or 
saturated (when available) thermal pressure. Plots of Log (akin) , Log (a mag ) , and 
Log(a to t) v s Log((3) are well fit by straight lines even as a^n, a ma9 ,and a to t vary 
by four orders of magnitude over the simulations included. The ratio akin/(y. mag 
and the product a to tl3 are quite constant and largely independent of the presence 
or absence of weak mean fields, the choice of initial and boundary conditions, and 
the resolution. In short, simulations have more strongly constrained the product 
a tot (5 than a to t itself. 
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1. Introduction 

Accretion disks are widely appreciated to be a source of emission from gas or plasma 
orbiting central stars or compact objects (c.f. Frank, King, Raine 1992). In order to explain 
the rapid variability and short lifetimes of accreting systems without unphysical mass densi- 
ties, some enhanced angular momentum transport beyond that which can be supplied by the 
microphysical transport coefficients is typically required. For sufficiently ionized disks, the 
magneto-rotational instability (MRI) offers a solution to this problem for sufficiently ionized 
disks (e.g. Balbus & Hawley 1991; 1998). 

The MRI feeds off of an initially weak magnetic field and the turbulence induced by the 
ensuing instability amplifies the fluctuating magnetic energy by line stretching. Sustained 
magnetic fields under the influence of a shear flow in a radially decreasing angular veloc- 
ity profile produce a negative magnetic (Maxwell) stress, which, in principle, produces the 
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dominant positive outward angular momentum transport. 3-D Simulations (e.g. Hawley, 
Gammie, Balbus, 1995,1996; Brandenburg et al. 1995; Stone et al. 1995) have revealed that 
the nonlinear evolution of systems unstable to the MRI leads to a Maxwell stress whose mag- 
nitude is larger than the negatively signed Reynolds stress. The MRI sustains the turbulent 
Maxwell stress and thus the outward angular momentum transport. 

While the MRI has been numerically shown to provide an effective turbulent magnetic 
stress, incorporating the saturated state of the MRI into the framework of practical accretion 
disc modeling using, for example the a to t viscosity coefficient formalism (Shakura & Sunyaev 
1973), where atot is defined from the turbulent viscosity, vt = 0tt o tC s H and c s and H are 
the sound speed and disk scale height) suffers from the non-universality of values of a to t 
inferred from simulations. Depending on the boundary conditions, initial conditions, different 
treatments of viscosity and resistivity, the presence or absence of stratification, and resolution 
(Pessah et al. 2007), the values inferred from simulation can vary by 4 orders of magnitude 
(see Figs. 1 and 2 below). However, the a to t prescription provides a practical mean field 
formalism that allows a straightforward calculation of accretion disk spectra for comparison 
to observations by parameterizing nonlinear correlations of turbulent fluctuations by a simple 
closure. Developing an improved mean field theory that also incorporates the physics of the 
MRI, while still being practical is an important target of recent and ongoing work (e.g. 
Ogilvie 2003; Pessah 2006ab). 

Here we emphasize that the ratio of the thermal to magnetic pressure, j3 is not generally 
an independent function of a tot , even though it is sometimes assumed to be in phenomeno- 
logical analytic disc models (e.g. Yuan et al. 2005). In this paper we combine the published 
data extracted into Tables 1-10 to determine the empirical correlation between the kinetic 
and magnetic contributions to a to t- 

In Sec. 2 we derive the formalism that relates the kinetic and magnetic parts of a to t to 
(3 for different adiabatic indices and give a physical argument for an inverse relation between 
a to t and (3. We do not present a rigorous theory in the present work as our main focus is 
empirical. Toward this end, in Sec. 3 we plot the data points from published simulations 
and infer the empirical values for the quantities defined in Sec. 2. The data reveal that the 
product atotP is nearly constant. We conclude in Sec. 4. 
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2. Maxwell and Reynolds Contributions to Transport and Relations to f3 

In the steady-state, ignoring microphysical viscosity, the mean azimuthal momentum 
equation is given by (e.g. Balbus & Hawley 1998) 
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The quantity inside the brackets represents the flux of angular momentum. Of particular 
interest is the r0 component of this flux, which, when greater than zero, represents the 
outward radial transport of angular momentum. It is given by 



B 4> D 

rpv^v r -r-^-B r 



(2) 



Because axisymmetric accretion disk equations formally represent mean field equations, 
we are interested in the averaged value of F r ^. Toward obtaining this, we split the magnetic 
field and velocity into mean (indicated by an overbar) and fluctuating components (indicated 
by lower case). As in Balbus & Hawley (1998), we take the mean to represent a height 
integration over all z, an average over all and an average over some fixed range of r. For 
a quantity Q = Q + q we have (q) = and 

(ffl = g= J %*£* . (3) 

where £ = J^ H pdz. Assuming that p = p (no fluctuating density), applying to (j2J) gives 

Xj?" 

F r <t> = \y<t>Vr - V A ^V A ,r + {v<t>V r ) ~ (& A)] ' ( 4 ) 

where is the Alfven velocity associated with the mean field and b is the Alfven velocity 
associated with the fluctuating field. The first term on the right is an inward flux of angular 
momentum, since V r < and V $ > at the inner most radii for an accretion disc. The 
remaining terms must provide the needed outward transport of angular momentum if matter 
is to accrete. The last two terms represent purely turbulent transport. In what follows, we 
assume that the mean magnetic field of smaller magnitude than the fluctuating field in 
saturation and that the dominant angular momentum transport comes from the last two 
terms of ([!]). (This is consistent with all of the simulations we consider.) 

The shearing box simulations of Table 1 employ local Cartesian coordinates in the 
rotating frame. In this shearing-sheet approximation, the mean velocity V y vanishes at the 
inner most radius r of the shearing box, and points in the — y for x = r — r > 0, decreasing 
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outward in Keplerian fashion such that V y = (r — r )r|^ ~ — Here £1 is the local 

orbital speed, Q is the orbital speed of the rotating frame, x = r — r , and x << r . In this 
context, we can combine the last 2 contributions of (jlj) into a Cartesian stress tensor 

W xy = (v y v x ) - (b y b x ). (5) 

Using the Shakura-Sunayev prescription of = o> to tC s H of Sec. 1, the stress for a Keplerian 
flow in a shearing box that corresponds to an outward flux of angular momentum is 

3 

- Veffd x Vy ~ -VLa to tC s H. (6) 

Setting this equal to ([SD gives a closure for the stress tensor, so that 

2W xy _ 2f(T)W xy 
tot ~3{lcM~ 3c? ' [) 



where T is the polytropic index and /(l) = y | (isothermal) and /(5/3) = y | (adiabatic). 
The last relation in (|7j) comes from solving the equation of hydrostatic equilibrium, namely 

1 dP GM z 



pdz ~~ R 2 R 

Over a density scale height, the solution gives QH = y^y, with midplane sound speed c s . 
We now split (J7j) into magnetic and kinetic terms such that a tot = akin + OL mag where 

= 2/(r)(W _ c ma9 (r,/3) 
mas " 3 C 2 /3 ' iyj 

and 

a, m = ^ = p , (10) 

where C mag (T, (3) = —^ip- ^^pf and Cfci n (r, /3) = ^rgzy- are to be inferred from the data. 
We also define C tot (r, (5) such that a to t = C to t/P- The statistically determined C to t need not 
exactly equal the separately determined best fit values of C mag and Cfc in . 

We can provide a crude physical argument which anticipates a strong anti-correlation 
between a to t and (3: Note that u e ff = a to tC s H = vtL, where vt and L are a turbulent 
velocity and dominant energy containing fluctuation scale. In a turbulent flow, the ratio of 
magnetic to kinetic turbulent energies is typically of order unity in saturation (and actually 
slightly larger than unity for MRI simulations). Crudely, if va ~ vt, then a to tC s H ~ vaL. 
But L ~ vt/Q ~ va/£1, if the eddy turnover time scale is comparable to the orbit time. The 
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latter is a reasonable assumption since the growth rate for a MRI instability that initiates 
turbulence is of order the rotation rate. We therefore have a to tC s H ~ v\/Q which implies 
OL to t ~ ^t/P' using the relations below Eq. (jSj). The specific anisotropy due to Keplerian 
shear likely implies a missing factor of order unity. 

Note that a mag , au n &tot, and (3 were defined above as a function of the actual thermal 
pressure at the instant of measurement, not the initial thermal pressure of a system from 
which the system could evolve. For the isothermal case, the distinction is not important 
because the pressure is a constant in time in the simulations, but for the adiabatic case the 
pressure evolves. For most of the published simulations, the transport coefficients and /3 are 
expressed with respect to the initial pressure at the midplane P = P(t — 0, z — 0), where 
the subscript indicates both the initial time and the midplane (z — 0). We distinguish this 
from the saturated thermal pressure at the midplane -P(O) = P(t » 0, z = 0). For Sano 
& Stone (2002), data for of the transport coefficients and (3 were published. Subsequently, 
we will explicitly distinguish P and P(0) where necessary. For the isothermal case we note 
that P = P(0) = P(z = 0) for all time. 



3. The a{(3) Relation from Published Numerical Simulation Data 

We have extracted data from published shearing box simulations to produce Tables 1- 
10. Note that as emphasized in the table captions, the double brackets used therein indicate 
the combination of a spatial average (as employed above) along with a time average at late 
tmes in the simulations. 

To make the numerical coefficients from the published data correspond with those of 
the previous section for in a mag and a^in we have multiplied the values from those references 
by 2 "CjT in Tables 1-10. Note also that the convention used in Brandenburg 1995 (see 
Brandenburg 1998) is a factor of -^jjj^ ~ 0.48 smaller than the convention used in the other 
references. Therefore, to match our numerical coefficients we have multiplied the a from 
that reference by by ^^p-j^T) = /TiTr ^° construct our Table 1 for a mag and a.u n - Also, the 
value of E mag in Brandenburg et al. (1995) must be divided by a factor of 2 when placed 
into our definition of /3 for the magnetic pressure. 

The simulations of Brandenburg et al. (1995) also differ from the others in that vertical 
field boundaries were used at the top and bottom of the box. This allowed mean field 
amplification, in contrast to the periodic boundary conditions used in the other simulations. 
However, since the mean field saturates at values relatively small compared to the random 
field, the effect of mean field growth on the total stress should be relatively small. 
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In Fig. labc we plot respectively, Log otkin(P(0)), Log a mag (P(0)), and Log a tot (P(0)) 
vs. Log /3(P(0)) for the T — 1 isothermal simulations, and plot Log afc in (P ) vs - Log a mag (P(0)) 
in Fig. Id. In Fig. 2 we show the analogous plots for the adiabatic T = 5/3 cases. There 
the distinction between P(0) and P is necessary. In Fig 2 we include the Sano et al. (2002) 
data points using P , whereas in Fig. 3 we give the analogous plots using P rather than P , 
for the cases in which au n {P)i «ma g (P), and /3(P) values are available. Although the use 
of P vs. Pq does not change products of the form af3, the range of values over which the x 
and y axes range can differ. From Fig. 3 however, it can be seen that the best fit curves 
still show a tight correlation and that the x and y axes range over more than an order of 
magnitude supporting the same general conclusions as the previous two Figs. 

The best fit solid lines are shown, with the equations for these lines at the top of 
each panel in Figs.l, 2, and 3. For each panel, the fit at the top is of the form Log(a) = 
ALog((3) + D, shown at the top with A and D constants. We can use these fits to extract 
best fit values of C mag (T,P) and Cfcj n (r, (3) defined in Sec. 2. Counter clockwise from the top 
right panel of Fig. 1 we then have 

c mag (i) = o.i3/r a03 , (ii) 

C Hn {l) = 0.02/? - 12 (12) 



c tot (i) = o.n/r - 06 , (13) 



and for Fig 2. we have 



C ma ,(5/3)=0.10/3(P )- a03 , (14) 

C fcm (5/3) = 0.02/3(P )°- 01 (15) 

CU(5/3) = 0.08/3(P )- a07 , (16) 
where we have added the explicit functional dependence on P . 
From Fig. 3 we have 

c mag (5/3) = o.n/r - 03 , (17) 

C«„(5/3) = 0.01/? - 20 (18) 

and 

CU(5/3) = 0.14/r a23 . (19) 



The results for Fig 3 are based on only 8 points, yet still the weak dependence of the 
near constancy of C mag and C to t is revealed in fTiT|) and f|T9|) . More robust are f|T3|) and fTT6l) 
which highlight a tight correlation. The results imply that C toi (l) ~ C iot (5/3)(P ) ~ 0.1 
even as a to t{Po) varies by more than 4 orders of magnitude. 
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The bottom right panels in Figs. 1 and 2 also highlight strict correlations between a^n 
vs. a mag . Fig Id shows that the best fit for the isothermal data is a^n = 0.22a£^ and 
Fig. 2d shows that for the adiabatic case, akm(Po) — 0.28a mag (Po)- The near constancy of 
otkinl&mag was also noted by (Pessah et al. 2006ab). 



4. Conclusion 

We have highlighted that data from published shearing box MRI simulations (Tables 
1-10) show tight anti-correlations of (3 with a mag , akin, and their sum atot- In particular, 
using our definitions, the product a tot (3 ~ 0.1, even as a to t varies by over 4 orders of mag- 
nitude between simulations. The data were taken from simulations invoking different codes, 
different vertical boundary conditions, different initial conditions, the presence or absence of 
stratification, the presence or absence of initial mean fields, explicit vs. numerical viscosity, 
different polytropic indices, and different resolutions. The data reveal that simulations have 
much more strongly constrained the product a to tP than any particular value of a tot . 

The fact that no universal value of atot emerges implies that the boundary and initial 
conditions, as well as the resolution (Pessah et al. 2007) are influencing its value. Pessah 
et al. (2007) show the magnitude of the magnetic stress correlates strongly with the box 
size. This is troubling as it implies no robust value of the transport coefficient presently 
emerges from simulations. In contrast, the narrow range of values for the constant product 
atotfl then must likely emerge from features identical in all simulations, such as the Keplerian 
shear profile This is consistent with expectations from analysis of the linear regime (Pessah 
et al. 2006a), and more general nonlinear closures (e.g. Ogilvie 2003; Pessah et al. 2006b). 
More work is needed to determine the constant a to t(3 from first principles. Alternatively, the 
constants atot/3 and a ma g/dkin can be employed as a constraints for closures. 

Modelers appealing to the MRI and extracting guidance from simulations should treat 
a to t and j3 as dependent parameters. The value of a to tf3 extracted from Fig. 2 emerges as 
much more robust universal constraint, than any specific value of a tot . A relation between 
a to t and (5 has been incorporated into some some disk models (e.g. Narayan et al. 1998) but 
not others (e.g. Yuan et al. 2005). 

Note also that the data herein are for thin disks. While similar principles would apply 
for thick disks, the numerical constants could be different. 
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Table 1. Brandenburg et al. 1995 



Run 


((SvrPo/^ 2 )) 


(2f(T)/3T)((-B x B y /4nP )} 


(2/(r)/3r)((p^/P » 


A 


154 


0.000734847 


0.000244949 


B 


143 


0.000342929 


9.79796E-05 


C 


40 


0.002008582 


0.000489898 


D 


143 


0.000685857 


0.000146969 


E2 


118 


0.001126765 


0.000146969 


AD 


40 


0.001518684 


0.000342929 



Note. — Runs with adiabatic equations of state and vertical initial fields. The 
double brackets indicate that all data are time- and volume-averages at late times 
over the turbulent layer of the disk. The magnetic field energy, Maxwell stress, 
and Reynolds stress are all normalized with respect to the initial gas pressure 
at the midplane, Pq. The run AD includes ambipolar diffusion. 
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Table 2. Fleming et al. 2000 



Run 


((SttPo/S 2 )) 


(2f(T)/3T)((-B x B y /4nP )} 


(2/(r)/3r)((p*avp » 


BZ1 


2.2222E+00 


5.8659E-02 


1.2009E-02 


BZ2 


3.6101E+00 


3.9260E-02 


9.2376E-03 


BZ3 


7.1429E+00 


2.0092E-02 


5.0807E-03 


BZ4 


1.6129E+01 


9.6995E-03 


2.7713E-03 


ZN1 


1.0000E+02 


1.0392E-03 


5.3116E-04 


ZN2 


3.3333E+02 


2.3094E-04 


2.3094E-04 


ZN3 


3.3333E+02 


2.0785E-04 


1.5935E-04 


Yl 


1.4184E+01 


6.9282E-03 


1.7782E-03 


Y2 


1.6393E+01 


6.0044E-03 


8.3138E-04 


Y3 


2.0408E+01 


2.3556E-03 


1.5935E-04 



Note. - Runs with adiabatic equations of state. Vertical field runs are 
labeled by the prefix BZ, zero net z runs are labeled by the prefix ZN, and 
toroidal field runs are labeled by the prefix Y. The double brackets indicate that 
all data are time- and volume-averages at late times over the turbulent layer of 
the disk. The magnetic field energy, Maxwell stress, and Reynolds stress are all 
normalized with respect to the initial gas pressure at the midplane, Pq. These 
runs include the effects of Ohmic resistivity. 
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Table 3. Hawley, Gammie, Balbus 1995 



Run 


((MPnlB 2 )) 


f2ffn/3n((-S T R,/47rPn)) 
\ J \ ii J \ \ x y I ^ ' 1 - 1 u// 


(2f(T)/3T)((pvJu 


Z3 


0.0294 


2.678905249 


0.743396207 


Z4 


1.9608 


0.069974853 


0.016165808 


Z5 


2.6385 


0.034410076 


0.013625466 


Z7 


2.5641 


0.050575884 


0.014780167 


Z9 


0.9346 


0.121474497 


0.027019993 


Z12 


1.0811 


0.110851252 


0.030022214 


Z15 


0.7686 


0.154037052 


0.032331615 


Z17 


3.9370 


0.029329394 


0.006928203 


Z18 


2.4038 


0.061430069 


0.027019993 


Z19 


2.1505 


0.066741691 


0.018475209 


Z20 


7.6336 


0.014087347 


0.003925982 


Z21 


6.5359 


0.017089568 


0.004618802 


Z22 


2.5974 


0.049190243 


0.010623245 


Z23 


27.0270 


0.004618802 


0.00092376 


Z24 


76.9231 


0.001385641 


0.00023094 


Z25 


30.3030 


0.003233162 


0.001154701 


Yl 


13.8889 


0.007159143 


0.002309401 


Y2 


30.3030 


0.003002221 


0.00092376 


Y3 


8.6207 


0.011316065 


0.003695042 


Y6 


23.8095 


0.003002221 


0.001154701 


Y7 


9.2593 


0.008544784 


0.002309401 


Y8 


1.6260 


0.019398969 


0.005080682 


Y9 


8.6207 


0.009468544 


0.002309401 


Y10 


62.5000 


0.00069282 


0.00023094 


Yll 


16.9492 


0.005542563 


0.001847521 


Y12 


50.0000 


0.001847521 


0.00069282 


Y13 


14.0845 


0.006928203 


0.002078461 


Y15 


8.5470 


0.010854185 


0.003233162 


Y17 


185.185 


0.000508068 


9.2376E-05 


Y18 


3.8610 


0.022632131 


0.006235383 


YZ1 


2.9674 


0.019860849 


0.010392305 
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Table 3 — Continued 



Run ({8ttP /B 2 }} 


(2/(r)/3r)<<-B x fl y /47rPo)) 


(2f(T)/3T)((pv x 5v y /P )) 


YZ2 3.9370 


0.026789052 


0.006235383 



Note. - Runs with adiabatic equations of state. Vertical field runs are 
labeled by the prefix Z and toroidal field runs are labeled by the prefix Y. 
Combined toroidal and vertical field runs are labeled by the prefix YZ. The 
double brackets inidcate that all data are time- and volume-averages at late 
times over the turbulent layer of the disk. The magnetic field energy, Maxwell 
stress, and Reynolds stress are all normalized with respect to the initial gas 
pressure at the midplane, Pq. 



Table 4. Hawley, Gammie, and Balbus 1996 



Run 


((SvrPo/^ 2 )) 


(2f(T)/3T)((-B x B y /4nP )} 


(2/(r)/3r)((p^/P » 


Rl 


34 


0.002817469 


0.000969948 


R2 


63 


0.001524205 


0.000623538 


R3 


32 


0.003140785 


0.000969948 


R4 


67 


0.001408735 


0.00057735 


R6 


16 


0.006189195 


0.002332495 


R7 


250 


0.000392598 


0.000161658 



Note. — Runs with adiabatic equations of state and random initial B. The 
double brackets indicate that all data are time- and volume-averages at late 
times over the turbulent layer of the disk. The magnetic field energy, Maxwell 
stress, and Reynolds stress are all normalized with respect to the initial gas 
pressure at the midplane, Pq. 
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Table 5. Sano k Stone 2002 



Run 


((SvrPo/5 2 )) 


/ / -B x B y 2/(r) v \ 
\\ 4nP 3r // 


({pv x 5v y /P )) 2 -§± 


Z02 


2.8094E+00 


3.7412E-02 


7.8058E-03 


Z03 


4.2117E+00 


2.4942E-02 


5.0345E-03 


Z04 


1.1940E+01 


9.0990E-03 


2.8637E-03 


SOI 


2.1511E+01 


3.8798E-03 


3.0022E-03 


S03 


7.7320E+00 


1.2009E-02 


3.0715E-03 


S10 


1.1988E+01 


7.5286E-03 


2.0346E-03 


Y02 


1.2686E+01 


8.1522E-03 


1.9237E-03 


Y04 


7.9931E+00 


1.2494E-02 


2.8868E-03 



Note. — Runs with adiabatic equations of state. Vertical 
field runs are labeled by the prefix Z, zero net z field runs are 
labeled by the prefix S, and vertical field runs are labeled by 
the prefix Y. The double brackets indicate that all data are 
time- and volume-averages at late times over the turbulent 
layer of the disk. The magnetic field energy, Maxwell stress, 
and Reynolds stress are all normalized with respect to the 
initial gas pressure at the midplane, P . These runs include 
the effects of the Hall term and Ohmic dissipation. 



- 15 - 



Table 6. Sano & Stone 2002 - saturated pressure 



Run 


((P/Po)) 


((P))/((By8n)) 


((-B x B y ))/((P))- 


((pv x 5v y ))/((P)) ■ 2 -§l 


Z02 


2.7800E+00 


7.8100E+00 


1.3458E-02 


2.8078E-03 


Z03 


4.6300E+00 


1.9500E+01 


5.3869E-03 


1.0874E-03 


Z04 


2.6800E+00 


3.2000E+01 


3.3952E-03 


1.0685E-03 


SOI 


2.5800E+00 


5.5500E+01 


1.5038E-03 


1.1637E-03 


S03 


1.9400E+00 


1.5000E+01 


6.1901E-03 


1.5832E-03 


S10 


3.3200E+00 


3.9800E+01 


2.2677E-03 


6.1283E-04 


Y02 


6.4400E+00 


8.1700E+01 


1.2659E-03 


2.9872E-04 


Y04 


5.8300E+00 


4.6600E+01 


2.1430E-03 


4.9515E-04 



Note. — Runs with adiabatic equations of state. Vertical field runs are labeled by 
the prefix Z, zero net z field runs are labeled by the prefix S, and vertical field runs 
are labeled by the prefix Y. The double brackets indicate that all data are time- and 
volume-averages at late times over the turbulent layer of the disk. In this table, P 
denotes the actual pressure at the time and position of averaging. These runs include 
the effects of the Hall term and Ohmic dissipation. 
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Table 7. Stone et al. 1996 



Run 


((SvrPo/5 2 )) 


(2/(r)/3r)<<^>> 


(2/(r)/3r)((p^/P )) 


AZ1 


42.2 


0.001769001 


0.000397217 


AZ6 


27.2 


0.00327935 


0.000739008 


AY1 


21.8 


0.002817469 


0.000672036 


AL6 


58.5 


0.001129297 


0.000274819 



Note. — Runs with adiabatic equations of state. Zero net z field runs 
are labeled by the prefix Z. AY1 is a vertical field run and AL6 has a 
"flux loops" initial field configuration. The double brackets indicate that 
all data are time- and volume-averages at late times over the turbulent 
layer of the disk. The magnetic field energy, Maxwell stress, and Reynolds 
stress are all normalized with respect to the initial gas pressure at the 
midplane, P . 
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Table 8. Fleming & Stone 2003 



Run 


«8ttP(0)/5 2 » 


(2f(T)/3T)((-B x B y /AnP(0))) 


(2/(r)/3r)(( P ^/p(o))) 


Zl 


6.6667E+01 


2.3570E-03 


4.7140E-04 


Z2 


3.3333E+02 


4.7140E-04 


1.4142E-04 


Z3 


1.0000E+03 


2.8284E-04 


4.7140E-05 


Z4 


3.5714E+01 


4.2426E-03 


1.4142E-03 


Yl 


5.0000E+02 


3.2998E-04 


4.7140E-05 


Y2 


8.3333E+01 


1.8856E-03 


4.7140E-04 



Note. — Runs with isothermal equations of state. Vertical field runs are labeled by 
the prefix Z and toroidal field runs are labeled by the prefix Y. The double brackets 
indicate that all data are time- and volume-averages at late times over the turbulent 
layer of the disk: 0.4 < \z/H\ < 2. (These models also include a central "dead" zone 
(0 < \z/H\ < 0.4) which is not included in the averaging.) The magnetic field energy, 
Maxwell stress, and Reynolds stress are all normalized with respect to the midplane 
gas pressure, -P(O). In these models, the ionization fraction varies with height. 
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Table 9. Miller & Stone 2000 



Run 


«8vrP(0)/5 2 » 


(2/(r)/3r)((- J B x S s /47rP(0)» 


(2f(T)/3T)((pv x 5v y /P(0))) 


BY1 


1.1547E+01 


1.1125E-02 


2.7483E-03 


DVO 
Dl A 


1.0 1 1 ori/+Ul 


o.oyiUri/-Uo 


A.VOOofj-Uo 


BY3 


3.8314E+01 


3.3564E-03 


1.0324E-03 


BY4 


1.7036E+01 


7.2596E-03 


2.0600E-03 


BY5 


1.8832E+01 


6.4582E-03 


1.6546E-03 


BY6 


1.0627E+01 


1.1597E-02 


2.8426E-03 


BYR1 


1.7575E+01 


6.7411E-03 


1.8243E-03 


BYR2 


1.8018E+01 


5.2326E-03 


1.3435E-03 


BYR3 


9.1743E+00 


1.6641E-02 


4.0447E-03 


ZN1 


4.5455E+01 


2.3099E-04 


1.0842E-04 


ZN2 


2.6882E+02 


2.5126E-04 


1.2681E-04 


BZ1 


5.5804E-01 


4.2756E-01 


2.6823E-02 


BZ2 


1.1351E+00 


1.4519E-01 


1.6735E-02 


BZ3 


9.6712E-01 


1.3114E-01 


1.4802E-02 



Note. - - Runs with isothermal equations of state. Toroidal field runs are labeled by 
the prefix BY, zero net z field runs are labeled by the prefix ZN, and pure z field runs are 
labeled by the prefix BZ. "R" identifies the resistive runs. The double brackets indicate 
that all data are time- and volume-averages at late times over the turbulent layer of 
the disk: < \z/H\ < 2. (These models also include an extended, quiescent corona 
(2 < \Z/H\ < 4) which is not included in the averaging.) The magnetic field energy, 
Maxwell stress, and Reynolds stress are all normalized with respect to the midplane gas 
pressure, -P(O). These runs include the effects of Maxwell's displacement current. 
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Table 10. Stone et al. 1996 



Run 


«87rP(0)/B 2 » 


(2/(r)/3r)«-B x V4^(o)» 


(2/(r)/3r)((p^/P(0))) 


IZl 


6.54E+01 


2.0930E-03 


5.8926E-04 


IZ2 


2.50E+02 


4.5726E-04 


2.1213E-04 


IZ3 


6.29E+01 


2.0836E-03 


5.9397E-04 


IZ6 


4.98E+01 


3.2150E-03 


8.9567E-04 


IY1 


2.79E+01 


4.6716E-03 


1.1314E-03 


IY2 


1.63E+02 


4.4312E-04 


1.4991E-04 


IY3 


2.58E+01 


3.7618E-03 


9.7109E-04 



Note. — Runs with isothermal equations of state. Toroidal field runs are labeled by 
the prefix IY and zero net z field runs are labeled by the prefix IZ. The double brackets 
indicate that all data are time- and volume-averages over the disk at late times. The 
magnetic field energy, Maxwell stress, and Reynolds stress are all normalized with 
respect to the midplane gas pressure, -P(O). 
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0.1 1 10 100 1000 10000 1e-04 0.001 0.01 0.1 1 

«8n P(0)/B 2 » (2f(r)/3r) «-B x B y /47iP(0)» 

Fig. 1. — Isothermal data of Tables 8-10. and best fit lines. Top row: Log afej n (P(0)) and 
Log a mag (P(0)) vs. Log /3(P(0)) respectively. Bottom row: Log a tat (P(0)) vs. Log /3(P(0)) 
and Log a)-i n (P(0)) vs. Log a mag (P(0)). The double brackets indicate that the data rep- 
resent a combination of spatial average and late time average (e.g. after 15 orbits). The 
symbols indicate the following specific data sets respectively: * = Fleming and Stone (2003); 
□ = Miller and Stone (2000); + = Stone et al. (1996). Values of the magnetic field energy, 
Maxwell stress, and Reynolds stress are all normalized with respect to the midplane gas 
pressure, -P(O), and for all of these runs, P(0) = 5 ■ 10~ 7 . For the most part, despite the 
different initial and boundary conditions and wide ranges of a mag , akin, and a tot , products 
of form a(3 lie close to the best fit lines shown. The line equations are at the top of each 
panel. 
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0.01 0.1 1 10 100 1000 " 0.0001 0.001 0.01 0.1 1 10 

«8jc P /B 2 » (2f(r)/3r) «-B x B y /47iP » 

Fig. 2. — Adiabatic runs from Tables 1-5,7 and the best-fit lines. Top row: Log otki n {Po) 
and Log a mag (Po) vs Log (5{Pq) respectively. Bottom row: Log a to t(Po) and Log akin(Po) 
vs. Log a ma g{Po) The double bracket averages are as in Fig. 1. P$ is the initial gas pressure 
at the midplane, which differs from the final midplane pressure P(0) for the adiabatic case 
(unlike the isothermal case of Fig 1. where the two pressures are interchangable.). The 
symbols indicate the following data sets respectively: A = Brandenburg et al. (1995); V = 
Fleming, Stone, and Hawley (2000); = Hawley, Gammie, and Balbus (1995); x = Hawley, 
Gammie, and Balbus (1996); = Sano and Stone (2002) (includes a Hall term); + = Stone 
et al. (1996). Equations for the best fit lines are given. 
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Fig. 3. — Same quantities plotted as in Fig. 2 but for the available adiabatic runs in which 
the saturated midplane pressure P(0) is used rather than the initial midplane pressure Pq of 
the previous figure. All data here are from Sano & Stone (2002). 



